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Abstract 

In this article, we study the structure of the p meson in the framework 
of the coupled rainbow Schwinger-Dyson equation and ladder Bethe-Salpeter 
equation with a confining effective potential. The u and d quark propagators 
get significantly modified, the mass poles are absent in the timelike region, 
which implements confinement naturally. The Bethe-Salpeter amplitudes of 
the p meson center around zero momentum and extend to the energy scale 
I/"") , about q 2 = lGeV 2 , which happens to be the energy scale of chiral symmetry 

breaking, strong interactions in the infrared region result in bound state. 
■ The numerical results of the mass and decay constant of the p meson are in 

agreement with the experimental data. 
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Low energy nonperturbative properties of quantum chromodynamics (QCD) put 
forward a great challenge to physicists as the SU(3) gauge coupling at low energy 
scale does not allow perturbative calculations. Among the existing theoretical ap- 
proaches, the coupled rainbow Schwinger-Dyson equation (SDE) and ladder Bethe- 
Salpeter equation (BSE) are typical. They have been very successful in describing 
the long distance properties of low energy QCD and the QCD vacuum (for reviews, 
one can consult Refs.[U [21 El HI IS])- The SDE can naturally embody dynamical 
chiral symmetry breaking and confinement, which are two crucial features of QCD, 
although they correspond to two very different energy scales [51 [6]. The possible 
interplay of two dynamics lies in the following two facts: the first one is that at 
high energy scale the chiral massive quark dresses itself with a gluon cloud and 
quark- ant iquark pairs, and creates a constituent quark mass; the second one is the 
double role of the light pseudoscalar mesons, as both Nambu-Goldstone bosons and 
qq bound states. The BSE is a conventional approach in dealing with two-body 
relativistic bound state problems [TJ. From solutions of the BSE, we can obtain use- 
ful information about the bound state structure of the mesons and obtain powerful 
tests of the quark theory. 
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Many analytical and numerical calculations indicate that the coupled rainbow 
SDE and ladder BSE with phenomenological potentials can give model independent 
results and satisfactory values [U [21 El HI E] • The pseudoscalar mesons, especially the 
7r and K, have been studied extensively to understand the bound state structures 
from nonperturbative QCD (for example, Ref. [8], for more literature, one can consult 
Refs. [H [21 El El E] ) • Chiral symmetry and axial Ward identity play an important role, 
in chiral limit, the dominating Bethe-Salpeter amplitude (BSA) of the pseudoscalar 
mesons with the structure 75 is given by B(q 2 ) / fp, where B is the scalar part of the 
quark self-energy and fp is the decay constant of the pseudoscalar meson. We can 
perform many phenomenological analysis without solving the BSE explicitly. 

The situation is much worse for the vector mesons. There is a conserved vector 
current and the corresponding Ward identity relates the longitudinal (not transverse) 
part of the vector vertex with the quark propagator. We cannot obtain dynamical 
insight into the relevant degrees of freedom without solving the BSE directly. Fur- 
thermore, it is more difficult to solve the vector BSE than the pseudoscalar BSE 
due to a larger number of covariants and higher masses, which require extrapolation 
into larger domain of the complex q 2 plane [9]. 

The often used effective potential models are confining Dirac 5 function potential, 
Gaussian distribution potential and flat bottom potential (FBP) [TUl [Tlj IT2] . The 
FBP is a sum of Yukawa potentials, which not only satisfies gauge invariance, chiral 
invariance and fully relativistic covariance, but also suppresses the singular point 
which the Yukawa potential has. It displays dynamical chiral symmetry breaking and 
confinement [T3], and works well in understanding the QCD vacuum (for example, 
the quark condensate, mixed condensate and vacuum susceptibility) as well as the 
meson structures (for example, the electromagnetic form factors, radius and decay 

constants) P21 IHIISl- 
In this article, we take the point of view that the p meson is the qq bound state, 
and study its structure and decay constant in the framework of the coupled rainbow 
SDE and ladder BSE with the infrared modified FBP, which take advantages of 
both the Gaussian distribution potential and the FBP. 

The article is arranged as the followings: we introduce the infrared modified 
FBP in section II; in section III, IV and V, we solve the coupled rainbow SDE and 
ladder BSE, study analytical property of the quark propagators, finally obtain the 
mass and decay constant of the p meson; section VI is reserved for conclusion. 

2 Infrared Modified Flat Bottom Potential 

The present techniques in QCD calculation cannot give satisfactory analytical results 
for the infrared behavior of the two-point Green's function of the gluon. Infrared 
enhanced effective potential models have been phenomenologically quite successful 

One can use a gaussian distribution function to represent the infrared behavior 
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of the two-point Green's function of the gluorj^], 



GM) = , (1) 

which determines the quark- ant iquark interaction through a strength parameter 
zu and a range parameter A [TTj . This form is inspired by the Dirac 5 function 
potential (in other words the infrared dominated potential) used in Ref.[10j. which 
it approaches in limit A — > 0. The integral J d 4 qq 2n Gi(q 2 ) is finite for spacelike 
squared momentum q 2 . Such an infrared behavior can result in large dressing for the 
quark's Schwinger- Dyson functions (SDFs) A(q 2 ) and B(q 2 ) near q 2 = 0, the curves 
at the infrared region may be steep enough to forbid extrapolation to deep timelike 
region. When we introduce an extra factor, q 2 /A, the modified gaussian distribution 
can result in more flat curve near zero momentum. Furthermore, systematic studies 
with the coupled SDEs of the quark, gluon and ghost indicate this type of behavior at 
about q 2 = [IB] . We use the following modified gaussian distribution to represent 
the infrared behavior of the two-point Green's function of the gluon [TT] . 

w 2 _l_ w 2 q 2 _l_ 

A^ 6 A A^ A 6 A • (2) 
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In numerical calculation, the range parameter A is taken to be a/A w 0.62Ge^, 

^2 

the Gaussian type of function q 2 e~~% centers around q = 0.6GeV, and extends to 
about q = 1.2GeV. Systematic studies with the coupled SDEs indicate that the 
nonperturbative gluon propagator is greatly enhanced at about q = lGeV [IS]. The 
value a/A pa O.Q2GeV is more reasonable than the value a/A = 0.3GeV taken in 
Ref.[9]. 

For intermediate momentum, we take the FBP as the best approximation and 
neglect the contribution from perturbative QCD calculations as the strong coupling 
constant at high energy scale is very small. The FBP is a sum of Yukawa potentials 
which is an analogy to exchange of a series of particles and ghosts with different 

masses, 

n 

q 2 + (N + jr]) 2 

where iV stands for the minimum value of the masses, r\ is their mass difference, 
and dj is their relative coupling constant. Definition of momentum regions between 
infrared and intermediate momentum is about Aqcd — 200MeV, which is set up 
naturally by the minimum value of the masses iV = IAqcd- Certainly, there are 
some overlaps between those regions, in this way, we can guarantee continuity for the 
momentum. The FBP at energy N+jp with j > 3 extends to the perturbative region 



2 In this article, we use the metric 5^ = (1,1,1,1), {7^7^ + 7^7^} = 26^, the coordinate 
= (it, x), the momentum = (iE, 
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and exhibits some perturbative characters. The infrared modified FBP is supposed 
to embody a great deal of physical information about all momentum regions. 

Due to the particular condition we take for the FBP, there is no divergence in 
solving the SDE. In its three dimensional form, the FBP takes the following form, 



e -(N+jv)r 



V(r) = (4) 
In order to suppress the singular point at r = 0, we take the following conditions, 



V(0) = constant 
dV(0) d 2 V(0) d n V{0) 



dr dr 2 dr T 



0. (5) 



The cij can be determined by solving equations inferred from the flat bottom condi- 
tion. As in previous literature [12], [131 HU EH] > n is se t to be 9. The gluon propagator 
can be approximated by 

G{q 2 ) = G 1 (q 2 ) + G 2 (q 2 ) . (6) 



3 Schwinger-Dyson equation 

The SDE, in effect the functional Euler-Lagrange equation of quantum field the- 
ory, provides a natural framework for studying the nonperturbative properties of 
quark and gluon Green's functions. By studying the evolution behavior and ana- 
lytic structure of the dressed quark propagator, one can obtain valuable information 
about dynamical chiral symmetry breaking phenomenon and confinement. In the 
rainbow approximation, the SDE takes the following form, 

/d^k X a X a 
-j^yl^S{k) lu —G^{k-p), (7) 

where 

S-\p) = iA(p 2 )-f -p + B{p 2 ) = A(p 2 )[ij -p + mip 2 )} , (8) 
G, v {k) = (8^-^)G(k 2 ). (9) 

In this article, we assume that a Wick rotation to Euclidean variables is allowed, 
and rotate p, k into the Euclidean region analytically. Alternatively, one can derive 
the SDE from Euclidean path-integral formulation of the theory, and avoid possible 
difficulties in performing the Wick rotation [17]. The analytical structures of quark 
propagators have interesting information about confinement, we will revisit this 
subject in the section V. 
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4 Bethe-Salpeter equation 



The BSE is a conventional approach in dealing with two-body relativistic bound 
state problems [TJ. Precise knowledge about the quark structure of the p meson can 
result in better understanding of its property. In the following, we write down the 
ladder BSE of the p meson, 

S-\q + g )x(q, P)S-\q - - ) = ~ J tJ^X^ P)l»G» v {q - k) , (10) 

where P^ is four-momentum of the center of mass of the p meson, is relative four- 
momentum between the two quarks, 7^ is the bare quark-gluon vertex, and x{Qi P) 
is the BSA of the p meson. 

We can perform the Wick rotation analytically and continue q, k into Euclidean 
region H. The Euclidean BSA of the p meson can be decomposed as 

x(q,P) = e, x ,(q,P), (11) 
+ {q,-^j^}{P^ + 2 i PF 3 } 

+ {q,-^j^\{FA + iPF,-i4F 6 +^F 7 } , (12) 

due to Lorentz covariance. Here e M is the polarization vector of the p meson [18]. The 

1 

BSAs Fi(q, q- P,P) can be expanded in terms of Tchebychev polynomials T n 2 (cos 9), 

00 

Fi(q,q ■ P,P) = Y,* n Fnq,P)q n P n Thcos9) , (13) 


where 9 is the included angle between g M and P^. It is impossible to solve an infinite 
series of coupled equations of the F", we have to make truncation in one or the other 
ways. Numerical calculations indicate that taking only some terms with n = 0,1,2 



3 To avoid any possible difficulties in performing the Wick rotation, we can derive both the SDE 
and BSE from Euclidean path-integral formulation of the theory directly, then continue the four- 
momentum of the center of mass of the p meson into Minkowski spacetime analytically P 2 = —m 2 p . 
As far as only the numerical values are concerned, the two approaches are equal. 
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can give satisfactory results @, 



+ - ^j^} K + ip v PF i - i + \p,m } ■ (14) 

In solving the BSE, it is important to translate the BSAs F™ into the same 
dimension of mass to facilitate the calculation, 

F ° - A°F °, F 2 - A 2 F 2 , Ff - A 1 ^ , - A 3 F 2 \ 
F 3 ° - A 2 F 3 ° , F 4 ° - A 1 ^ , F5 - A 4 F 5 \ F 6 ° - A 2 F 6 °, 
F 7 ° -> A 3 F 7 °, g -> q/A, P -> P/A , 

where A is a quantity with dimension of mass. The ladder BSE of the p meson can 
be projected into the following nine coupled integral equations, 

"£H(i,j)F^ 2 (q,P) = ^ fd*kK(i,j), (15) 

where H(i,j) and K(i,j) are 9x9 matrices, the corresponding ones for the pseu- 
doscalar mesons are 4x4 matrices [H]. The analytical expressions of the matrix 
elements H(i,j) and K(i,j) are cumbersome and will take up more than seven 
pages, and not shown explicitly for simplicity. 

4 We can borrow some ideas from the twist- 2 light-cone distribution amplitudes <p((i, u) of the 
pseudoscalar mesons it and K, where u is momentum fraction of the quark, and \x is the energy 
scale. The <)>((i,u) is always expanded in terms of Gegenbauer polynomials, <fr(fj,,u) = 6u(l — 
u) {l + ai(M)C 1 3/2 (2u - 1) + a 2 (/i)C 2 3/2 (2w - 1) + a 4 (fi)Cl /2 (2u - 1) + • • ■ }, where C 3/2 (2u - 1), 

C^ 2 (2u — 1) and C^ 2 {2u — 1) are Gegenbauer polynomials, and aj(^t) are nonperturbative coef- 
ficients. The coefficients ai(ii) can be estimated with the QCD sum rules approach, for large i, 
the Oi(/x) involve high dimension vacuum condensates which are known poorly. In general, one 
can retain only the first few terms and fit them with experimental data, and the truncated light- 
cone distribution amplitude <p(fi, u) always gives satisfactory results |19) . In this article, we retain 

1/2 

only some terms with n = 0, 1 and 2 in expansion with Tchebychev polynomials T n (cos#), 
the contributions from other terms are supposed to be small and neglected here. If the contribu- 
tions from the neglected terms are large, they cannot decouple approximately from the BSEs and 
warrant the BSEs have reasonable solution. The numerical results indicate that the dominating 
contribution comes from the Fq, the sub-dominating contributions come from the F® and F®, 
and F§ ^> (frriiFQ, we expect the truncation is reasonable. It is obvious that if we take into 
account more terms in expansion, more accurate values can be obtained. However, the analytical 
expressions of the coupled BSEs become very clumsy and are beyond capability of our computer 
in numerical calculations. 
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We can introduce a parameter A(P 2 ) and solve above equations as an eigenvalue 
problem. If there really exist a bound state in the vector channel, the mass of the p 
meson can be determined by the condition X(P 2 = —m 2 p ) = 1, 



Y^H{i,j)Ff> lfl {q,P) = \(P 2 )J2 d 4 kK(i,j). (16) 
j j J 

The matrix elements H(i, j) are functions of the quark's SDFs A{q 2 + ^-±q- P) 
and B(q 2 + ^- ±q ■ P) . The relative four- momentum g M is a quantity in Euclidean 
spacetime, while the center of mass four-momentum P^ is a quantity in Minkowski 
spacetime, q ■ P varies throughout a complex domain. We can expand the A and B 
in terms of Taylor series of q ■ P to avoid solving the SDE with complex values of 
quark momentum, for example, 

A(q 2 + ^±q.P) = A{q 2 + ^)±A{q 2 + ^)'q.P + ... . 

The other problem is that we cannot solve the SDE in the timelike region. The two- 
point Green's function of the gluon cannot be exactly inferred from the SU(3) gauge 
theory even in small spacelike momentum region. We can extrapolate the values 
of the A(q 2 ) and B(q 2 ) from the spacelike region smoothly to the timelike region 

2 

q 2 = — j- pd — Q.lbGeV 2 with suitable polynomial. The masses of the vector mesons 
are larger than the pseudoscalar mesons. So it is very difficult to extrapolate the 
values to the deep timelike region. We must be careful in choosing the polynomial 
to avoid possible violation of confinement in sense of appearance of pole masses 
q 2 = —m 2 (q 2 ) in the timelike region [TTJ [TJ]. This requires a certain amount of 
fine tuning. Furthermore, if the curves of the A(q 2 ) and B(q 2 ) are very steep near 
q 2 = 0, very large values of the A(q 2 + ^-) and B(q 2 + ^-) are obtained in the 
timelike region {q 2 + ^- -< 0), the solution of the BSE in the small spacelike region 

(0 z< q 2 -< is not reasonable. In this article, we use the modified Gaussian 
distribution rather than the Gaussian distribution to modify the infrared behavior 
of the two-point green's function of the gluon to outcome the above difficulties. 
Finally we write down normalization condition for the BSAs of the p meson, 

f / p)s:1 + * s;lxfe p)d M] = 2p » • (i7) 

where x = 74X + 74, 5+ = S(q + ~) and 5_ = S(q — ~). 

5 Coupled rainbow SDE and ladder BSE 

Now we study the coupled equations of the rainbow SDE and ladder BSE of the p 
meson. 
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In order to demonstrate confinement of quarks, we take the Fourier transform 
with respect to Euclidean time T for the scalar part of the quark propagator [UGH [20], 



S* S (T) = -^e 



-oo 



dq* iQ4 T B(q ) 

2vr q 2 A 2 {q 2 ) + B 2 {q 2 ) lq=0 ' 1 ' 

where the 3-vector part of g M is set to zero. If the S(q) has a mass pole at q 2 = 
—m 2 (q 2 ) in the real timelike region, the Fourier transformed S*(T) would fall off 
as e~ mT for large T. In numerical calculations, for small T, the values of S* are 
positive, and decrease rapidly to zero with increase of T, which are compatible with 
the result (tendency of curve with respect to T) from lattice simulations [21]. For 
large T, the values of S* are negative, except for occasionally a very small fraction 
of positive values. The negative values of S* indicate an explicit violation of axiom 
of reflection positivity [22], the quarks are not physical observable. 

The u and d quarks have small current masses, the dressing or renormalization 
is large and the curves of the SDFs are steep, which indicates dynamical chiral 
symmetry breaking phenomenon. At zero momentum, m u (0) = rrid(0) = 0.51GeV, 
the Euclidean constituent quark masses are m u {m u ) = m^m^) = 0A2GeV, which 
are compatible with the constituent quark masses in literature. From the solutions 
of BSEs of the p meson as an eigenvalue problem, we obtain the mass, 

m p = 770MeV . (19) 

It is obvious 

m u {m u ) + m d (m d ) > m p . (20) 

The attractive interaction between the quark and antiquark in the infrared region 
can result in bound state. 

The dominating and sub-dominating components of the BSAs are the F®, F® 
and F®, the six components F 2 , F£, F®, F^, F® and F® are very small and of 
minor importance. From plotted BSAs (see Fig.l and Fig.2 as examples), we can 
see that the BSAs of the p meson have analogical momentum dependence, while 
the quantitative values are different from each other. Just like the qq , qQ and QQ 
pseudoscalar mesons [H], the BSAs of the p meson center around zero momentum 
and extend to the energy scale about q 2 = lGeV 2 , which happens to be the energy 
scale of chiral symmetry breaking. The strong interactions in the infrared region 
result in bound state. The BSAs of the p meson can give satisfactory value for the 
decay constant, which is defined by 

fpmpe^ = (0|g7 M g|p(P)) , 

= N e jTr[ W (k,P)]-^. (21) 

Carrying out trace explicitly, we can see that only the BSAs F® and F$ are relevant 
to the decay constant. The F$ is numerically very small (see Fig.2), the dominating 
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Figure 1: The dominating and sub- dominating components of the BSAs . 
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Figure 2: One of the small components of the BSAs 
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contribution comes from the Fq. Finally we obtain the value of the decay constant, 

f p = 223MeV , (22) 

which is in agreement with the experimental data. 

In calculation, the input parameters are taken as N = 1.0A, V(0) = —14. OA, 
7] = 5.0A, m u = m d = 6MeV, A = 200MeV, w = l.SGeV and A = 0.39GeV 2 . 

6 Conclusion 

In this article, we study structure of the p meson in the framework of the coupled 
rainbow SDE and ladder BSE with the confining effective potential. By solving the 
coupled rainbow SDE and ladder BSE as an eigenvalue problem numerically, we 
obtain the SDFs, BSAs, mass and decay constant of the p meson. 

The dressing (or renormalization) for the SDFs of the u and d quarks is large 
and the curves are steep, which indicate dynamical chiral symmetry breaking phe- 
nomenon. The mass poles are absent in the timelike region, which implements 
confinement naturally. The BSAs of the p meson have analogical momentum de- 
pendence, while the quantitative values are different from each other, where the 
dominating and sub-dominating components are the Fq, F° and F®, other six com- 
ponents Fq, F%, F®, F§, Fq, F® are of minor importance. The BSAs center around 
zero momentum and extend to the energy scale about q 2 = lGeV 2 , which hap- 
pens to be the energy scale of chiral symmetry breaking. Strong interactions in the 
infrared region result in bound state. 

The numerical results of the mass and decay constant of the p meson are in 
agreement with the experimental data. Once satisfactory SDFs and BSAs of the p 
meson are obtained, we can use them in other phenomenological analysis. 
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